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I. INTRODUCTION 


In nonleptonic hyperon decays such as A — p7r~, it is possible to search for CP violation by 
comparing the angular distribution with that of the corresponding anti-hyperon decay [1]. The 
Fermilab experiment HyperCP is currently analyzing data searching for CP violation in such a 
decay. 

The reaction of interest for HyperCP is the decay of a polarized A, with known polarization w, 
into a proton (whose polarization is not measured) and a tt” with momentum q. The interesting 
observable is a correlation in the decay distribution of the form 


dP 

dU 


1 + a w ■ q . 


( 1 ) 


The branching ratio for this mode is 63.9%, and the parameter a has been measured to be 
= 0.642 [2]. The CP violation in question involves a comparison of the parameter a with 
the corresponding parameter a from the reaction A — p7i~^. 

To obtain polarized A’s with known polarization, it is necessary to study the double decay chain 
Att” —> p7i~7r~ [3, 4]. This eventually leads to the experimental observable being sensitive 
to the sum of CP violation in the S decay and CP violation in the A decay. 

In both reactions, Att” and A — p7i~, the hnal state can be reached from the initial state 

via |A/| = i or |A/| = | transitions. It is known that due to the existence of a strong |A/| = | 
rule for nonleptonic hyperon decay, the dominant contribution to the CP-violating asymmetries 
arises from interference between an A-wave and a P-wave within the |A/| = | transition [5-7]. 
One can dehne the CP-violating asymmetries 


- ^(A°) 




As 


A(s;) 




for the A and S decays, respectively. The experimental observable is then [3, 4] 


( 2 ) 


Asa — A^ A= . (3) 

Approximate expressions have been obtained for in the case of |AJ| = ^ dominance [6], 
namely 

Aa ^ - tan (4 - sin (0^ - ) , 

„ „ „ „ (^) 

Ac ~ — tan {5p — (5|) sin — 0^) . 

Here, 5g [5p) is the strong P-wave (P-wave) Ntt scattering phase-shift at y/s = Af^, and (^5p) is 
the strong P-wave (P-wave) Att scattering phase-shift at y/s = Me. Moreover, (0p“) are the 
CP-violating weak phases induced by the |AP| = 1, |A/| = | interaction in the P-wave (P-wave) 
of the A — p'K~ and Att” decays, respectively. 

Experimentally, the current published limit is Ac^ = 0.012 ± 0.014 from E756 [3], and the 
expected sensitivity of HyperCP is 10“^ [4]. In addition, HyperCP has recently obtained a pre¬ 
liminary measurement of Ac^ = (—7 ± 12 ± 6.2) x 10“^ [8]. Previous estimates for Ac^ indicated 
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that it occurs at the few times 10“^ level within the standard model [6, 7, 9] and that it can be as 
large as 10“^ beyond the standard model [6, 10]. The larger asymmetries occur in models with an 
enhanced gluon-dipole operator that is parity-even and thus does not contribute to the e' parameter 
in kaon decay. The 10“^ upper bound corresponds to the phenomenological constraint from new 
contributions to the e parameter in kaon mixing. This illustrates the relevance of the HyperCP 
measurement which complements the e' experiments in the study of CP violation in IAS"! = 1 
transitions. 

The strong nN scattering phases needed in Eq. (4) have been measured to be Sg ~ 6° and 
6p ~ —1° with errors of about 1° [11]. In contrast, the strong An scattering phases have not been 
measured. Modern calculations based on chiral perturbation theory indicate that these phases are 
small, with |(5|| being at most 7° [12-17]. For our numerical estimates, we will allow the An phases 
to vary within the range obtained at next-to-leading order in heavy-baryon chiral perturbation 
theory [15], 

-3.0° < (Jf < +0.4° , -3.5° < < -1.2° . (5) 

One could choose to be less constrained and include the larger = —7° found in Ref. [15], but this 
would only enlarge the range and hence the uncertainty of the predicted asymmetry. In any case, 
eventually these phases can be extracted directly from the measurement of the decay distribution in 
S —An [4, 18, 19]. Recently E756 has reported a preliminary result of — = 3.17° ±5.45° [18]. 

In this paper, we estimate the weak phases that appear in and within the standard model. 
In Sec. II, we present a calculation of the weak phases guided by heavy-baryon chiral perturbation 
theory in terms of three unknown weak counterterms. In Sec. Ill, we estimate the value of these 
counterterms by considering contributions that arise from the factorization of the penguin operator 
and also nonfactorizable contributions estimated in the MIT bag model. Sec. IV contains the 
resulting weak phases and CP-violating asymmetries. Finally, in Sec. V, we compare our results 
to those of previous work and present our conclusions. For completeness, we also provide in an 
appendix the results for the corresponding asymmetries in S —Nn decays. 


II. CHIRAL PERTURBATION THEORY 


The chiral Lagrangian that describes the interactions of the lowest-lying mesons and baryons is 
written down in terms of the lightest meson-octet, baryon-octet, and baryon-decuplet fields [20-23]. 
The meson and baryon octets are collected into 3x3 matrices and P, respectively, and the 
decuplet helds are represented by the Rarita-Schwinger tensor which is completely symmetric 
in its SU(3) indices (a, b, c). The octet mesons enter through the exponential S = .^^ = exp(i(p//), 
where / is the pion-decay constant. 

In the heavy-baryon formalism [23, 24], the baryons in the chiral Lagrangian are described by 
velocity-dependent fields, and T^. For the strong interactions, the leading-order Lagrangian is 
given by [23-25] 

41 ) = + (B„ iv ■ CB„> + 2D {B„S; 2F {B,S; [.4„, B,] > 

- f"it. • + Am f+ C (f^A^B„ + B„An:) + 2H f"S. ■ AT,^ , ( 6 ) 
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where (■ ■ ■) denotes Tr(- ■ ■) in flavor-SU(3) space, is the spin operator, and 

d 

\ = Hi ) = ^ + 0(4’“) . (7) 

with further details given in Ref. [26]. In this Lagrangian, D, F, C, and are free parameters, 
which can be determined from hyperon semileptonic decays and from strong decays of the form 
T —> B(j). Fitting tree-level formulas, one extracts [23, 24] 

D = 0.80 , F = 0.50 , \C\ = 1.7 , (8) 

whereas 77 is undetermined from this £t. From the nonrelativistic quark models, one finds the 
relations [25] 


3F = 2D , C = -2D , n = -3D , (9) 

which are well satisfied by D, F, and C, suggesting the tree-level value 

n = -2.4 . (10) 

In our numerical estimates, we use Eqs. (8) and (10) for the leading-order results and the estimate 
of their uncertainty from one-loop contributions, with C and 77 only appearing in loop diagrams 
involving decuplet baryons. As another estimate of the uncertainty in these results, we will evaluate 
the effect of varying D and F between their tree-level values above and their one-loop values to be 
given later. 

At next-to-leading order, the strong Lagrangian contains a greater number of terms [27]. The 
ones of interest here are those that explicitly break chiral symmetry, containing one power of the 
quark-mass matrix M = diag(0, 0, . For our calculation of the factorization of the penguin 

operator, we will need these terms in the form 

= i/bx+) + ^(4{x+.S.}) + 4|-(4[x+.Bj) + ^(x+)(4B„) 

+ ^ ^ (x+> , (11) 

where we have used the notation x+ = to introduce coupling to external (pseudo)scalar 

sources, y = s -f fp, such that in the absence of the external sources y reduces to the mass matrix, 
X = 2BqM. As will be discussed in the next section, we also need from the meson sector the 
next-to-leading-order Lagrangian 

4“'= L5(a'‘Etsx«“x+?> + . (12) 

where only the relevant term is explicitly shown. In Eqs. (11) and (12), the constants Bq, b^pQ, c, 
Cq, and Lg are free parameters to be fixed from data. 

As is well known, the weak interactions responsible for hyperon nonleptonic decays are described 
by a IAS"! = 1 Hamiltonian that transforms as (8 l, 1r) © (27l, Ij^) under SU(3 )lXSU(3)j^ rotations. 
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It is also known from experiment that the octet term dominates the 27-plet term, as indicated by 
the fact that the |AJ| = | components of the decay amplitudes are larger than the |A/| = | 
components by about twenty times [26, 28]. We shall, therefore, assume in what follows that the 
decays are completely characterized by the (8 l, | AJ| = | interactions. The leading-order chiral 

Lagrangian for such interactions is [20, 29] 

+ lJ^hd^T.d‘‘T)) + -a.c. , (13) 

where h is a 3x3 matrix with elements parameters and 7 g contain the 

weak phases to be discussed below. 

The weak Lagrangian in Eq. (13) is thus the leading-order (in yPT) realization of the effective 
IAS"! = 1 Hamiltonian in the standard model [30], 

^ 10 

w. = ^ Kill, ^ C. Q, + H.c , (14) 


where Gp is the Fermi coupling constant, V],; are the elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [31], 


Ci = Zi + TVi = Zi 


VuVts 


(15) 


are the Wilson coefficients, and are four-quark operators whose expressions can be found in 
Ref. [30]. Later on, we will express V],; in the Wolfenstein parametrization [32]. It follows that 

VdV, = A , v:,v„ = -AM^ (1 - P + i,,) (16) 


at lowest order in A. For our numerical estimates, the relevant parameters that we will employ 
are [33] 


A = 0.2219 , A = 0.832 , r/ = 0.339 . 


( 17 ) 


In the next section, we match the penguin operator Qq in the short-distance Hamiltonian of Eq. (14) 
with the corresponding Lagrangian parameters in Eq. (13). 

We now have all the ingredients necessary to calculate the weak decay amplitudes in terms of the 
four parameters hp p^ and jg (only the hrst two are needed at leading order). In the heavy-baryon 
formalism, the amplitude for the weak decay of a spin-| baryon B into another spin-| baryon B' 
and a pseudoscalar meson 0 has the general form [29] 

= hg/ + 2S^-p^App,^'j Up , (18) 

where the superscripts refer to the S- and P-wave components of the amplitude. To express our 
results, we also adopt the notation [29] 

a/2 f 
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(19) 



u 




a 


(a) (b) 

FIG. 1: Leading-order diagrams for (a) 5-wave and (b) P-wave hyperon nonleptonic decays. In all figures, 
a solid (dashed) line denotes a baryon-octet (meson-octet) field, and a solid dot (hollow square) represents 
a strong (weak) vertex, with the strong vertices being generated by in Eq. (6). Here the weak vertices 
come from the hj^ p terms in Eq. (13). 


With the Lagrangians given above, one can derive the amplitudes at leading order, represented by 
the diagrams in Fig. 1. Fig. 1(a) indicates that the S'-wave is directly obtained from a weak vertex 
provided by Eq. (13). The leading contribution to the P-wave arises from baryon-pole diagrams, 
as in Fig. 1(b), which each involve a weak vertex from Eq. (13) and a strong vertex from Eq. (6). 
Thus the leading-order results for amplitudes not related by isospin are [20, 29] 

= -^{^D + ^hp) , = -^[hp-^hp) , 

(20a) 

- 0 - -h 

“s+n7r+ “ ’ “s-nTT- “ '''D + '^F i 


(P) 

'Ap 



a 


(p) 

E+n 7 r+ 


a 


(P) 

Yl~nn~ 


2D {hp — hp^ ^ (D -\- F) {hp -\- 3hp^ 

^/6 (m^ - m^) y/6 [m^ - m^) 

-2D {hp + hp) {D - F) [hp - 3hp) 

a/ 6 (m^ — rrij^ a/ 6 (m^ — 

_ -D {hp - hp) \D {hp + 3hp) 

_ -F [hp - hp) \D {hp + 3hp) 


(20b) 


The leading nonanalytic contributions to the amplitudes arise from one-loop diagrams, with hp 
only appearing in those involving decuplet baryons. These contributions have been calculated by 
various authors [20, 29, 34, 35], and we will adopt the results of Ref. [35] for the numerical estimate 
of our uncertainty. 

In Fig. 2, we show the kaon-pole diagram to be discussed later on. In this diagram, there is a 
strong vertex from Eq. (6) followed by a kaon pole and a weak vertex from the 7 g term in Eq. (13). 
Notice that this term is not only subleading in the chiral expansion, but also suppressed by an 
ni^/rn]^ factor (and hence vanishing in the = 0 limit). 

Once the values of the weak couplings hp p are specihed, the formulas in Eq. (20) determine the 
leading-order amplitudes. It is well known that this representation does not provide a good £t to the 
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FIG. 2; Kaon-pole diagram contributing to P-wave hyperon nonleptonic decays. The weak vertex here 
comes from the 7 g term in Eq. (13) 


measured P-wave amplitudes, and that higher-order terms are important [20, 22, 29, 34-36]. The 
procedure that we adopt for estimating the weak phases is to obtain the real part of the amplitudes 
from experiment (assuming no GP violation) and to use Eq. (20) to estimate the imaginary parts. 
The dominant GP-violating phases in the |AJ| = ^ sector of the jAPI = 1 weak interaction occur 
in the Wilson coefficient Gg associated with the penguin operator Qg. Our strategy will be to 
calculate within a model the imaginary part of the couplings hj^pQ and yg induced by Qq. As a 
numerical result, we propose a central value from leading-order yPT [Eq. (20)] and an estimate of 
the error from the nonanalytic corrections obtained with the expressions given in Ref. [35]. 


To end this section, for later use we collect in Table I the experimental values of the S- and 
P-wave amplitudes of interest, reproduced from Ref. [35]. The numbers are extracted (neglecting 
strong and weak phases) from the measured decay width T and decay parameter a by means of the 
relations 


T 


\Pb'\ 
4:71171 p 


[Epi + 


rriB') (|sP + bP) , 


2 Re(s*p) 
|sp -I- b|2 


( 21 ) 


The s and p amplitudes are related to those in Eq. (18) above by^ 


s = , p = \pp,\A^^'^ . 


( 22 ) 


TABLE I: Experimental values for S- and P-wave amplitudes, in units of Gpm^. 


Decay mode 

s 

P 

A- 

-^p7T~ 

1.42±0.01 

0.52±0.01 

- 


-1.98±0.01 

0.48±0.02 

E+- 


0.06±0.01 

1.81±0.01 

E+- 

->pTr^ 

-1.43±0.05 

1.17±0.06 

E"- 

-^n7T~ 

1.88±0.01 

-0.06±0.01 


^ In Refs. [29, 35], the p expression has the opposite sign, p = —\pp,\A^^\ but this turns out to be inconsistent 
with the amplitude formula from which both T and a are derived. Nevertheless, the sign flip does not affect 
the conclusions of Refs. [29, 35], as the fits therein were performed to the S'-waves and the P-waves were poorly 
reproduced regardless of the sign of p. 
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III. ESTIMATE OF COUNTERTERMS 


Our task in this section is to match the dominant | A/| = ^ UP-violating term from the standard- 
model effective weak Hamiltonian in Eq. (14) to the weak chiral Lagrangian in Eq. (13). That is, 
to compute the imaginary part of the parameters h^, hp, Hq, and 7 g that is induced by ImUgQg 
in Eq. (14). To do this, we will include both factorizable contributions, that arise from regard¬ 
ing the operator as the product of two (pseudo)scalar densities, and direct (nonfactorizable) 
contributions calculated in the MIT bag model. 

The nonfactorizable contributions are easily obtained from the observation that the weak chiral 
Lagrangian of Eq. (13) is responsible for nondiagonal “weak mass terms” such as 

{n\{n^)^\A) = u^uj^ , 

(A|(^w)8lS°> = , (23) 

I (^^w) 81 ^ ~ ^ ) 


where the subscript 8 denotes the component of that transforms as (8 l, Ij^). These terms can 
be computed directly from the short-distance Hamiltonian in Eq. (14) by calculating in the MIT 
bag model the baryon-baryon matrix elements of the four-quark operators. From the basic results 
in Appendix A, one hnds the contributions 


hp) — 


^ Ug (3 a — 5 6) , 


V2 




hp = Cf^ia+ 

v2 

Ug (—12 a-|-4 6) , 


(24) 


where a and b are bag parameters whose values are given in Eq. (All) for hp p and in Eq. (A12) 
for hp. Numerically, the imaginary part of Ug then yields, in units of f^Gpm1+ X^A^r], 


Imhp = 0.278 , Imh^ = 1.04 , Imhp = —3.13 |/g , (25) 


where ~ 92.4 MeV has been used. The units are chosen to separate both the conventional 
normalization for the hyperon decay amplitudes, as in Eq. (19) and Table I, and the relevant 
combination of CKM parameters occurring in the observables 

To obtain the factorizable contributions to the imaginary part of the parameters hp pp, we follow 
the procedure used in kaon physics for jg [37]. As shown in Appendix B, the lowest-order chiral 
realization of a factorized Qq contributes to the weak Lagrangian in Eq. (13) with 


hp - 

hp 


Up A 


8 Ug fB, bp , 


U,A 




V2 


8C,fB,bp, 


h 

T, = kU 16 C,BlL, . 


( 26 ) 



The values of bj^, bp, and c can be determined by fitting the mass formulas derived from the La- 
grangian in Eq. (11), with y = 2BqM, to the measured masses of the octet and decuplet baryons [2], 
Thus we find 

bprifig = 0.0301 GeV , bpirig = —0.0948 GeV , = 0.221 GeV (27) 

for = 0. In this limit, the Lagrangian in Eq. (11) also gives m\- = Using 

nig = rhg{fi = m^) = 170 MeV from Ref. [30], we then have 

bp = 0.177 , bp = -0.558 , c = 1.30 , Bq = 1.45 GeV . (28) 

Eor Lg, we adopt the value Lg = 1.4 x 10“^ found in Ref. [38]. Setting / = /^ ~ 92.4MeV, we 

then obtain the Qq contributions, in units of \/2 f^Gpm‘l+ rjX^, 


Imhp = 4.86 i/g , Imh^ = — 15.3|/g , Imh^ = 35.6 i/q , 

Im 78 = 18.8 , 


(29) 


where the formula with 7 g is the usual one appearing in the calculation of e' in kaon decay, and 
we have introduced the standard parameter to encode deviations from factorization [30], so 

that here = 1. 


IV. NUMERICAL RESULTS 

If Eq. (20) provided a good fit to the hyperon decay amplitudes, it would be straightforward to 
calculate the weak phases of Eq. (4). We would simply divide the imaginary part of the amplitudes 
by the real part of the amplitudes obtained from a matching of the parameters hp p to the short- 
distance Hamiltonian. However, as we mentioned before, leading-order chiral perturbation theory 
fails to reproduce simultaneously the S- and R-wave amplitudes. Gonsequently, we are forced to 
employ the real part of the amplitudes that are extracted from experiment under the assumption 
of no GR violation. 

An additional problem occurs if we calculate the imaginary part of the amplitudes from a match¬ 
ing of the full weak Hamiltonian to hp p and then divide it by the experimental amplitudes, as this 
introduces spurious phase differences. This can be easily understood by considering the case where 
only one operator occurs in the short-distance weak Hamiltonian. In such a case, it is clear that 
there can be no GR violation, as there is only one weak phase in the problem. However, if we use 
the procedure outlined above to calculate the phase difference (j)g — (j)p, we obtain a nonzero result 
due to the mismatch between the predicted and the measured ratio p/s. 

On the other hand, if there are two operators in the short-distance weak Hamiltonian, and 
one of them is mostly responsible for the real part of the amplitudes while the other one is mostly 
responsible for the weak phases, the procedure above does not introduce spurious phases. Of course, 
the predictions obtained are reliable only to the extent that the model reproduces the true imaginary 
part of the amplitudes. 
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In view of all this, we adopt the following prescription to obtain the weak phases. We hrst assume 
that the real part of the weak decay amplitudes originates predominantly in the tree-level operators 
Qig. This is true in the bag model, for example, as can be seen from the results in Appendix A. 
We then assume that the imaginary part of the amplitudes is primarily due to the ImC'gQg term 
in the weak Hamiltonian. This is true both in the bag model and in the vacuum-saturation model 
of Ref. [7], and is due to the purely |A/| = ^ nature of the CP observables With these 

assumptions, we calculate a central value for the imaginary part of the weak decay amplitudes 
using Eq. (20) with values for Im p obtained in the previous section by adding the factorizable 
and nonfactorizable contributions. We estimate the uncertainty in this prediction by computing 
the leading nonanalytic corrections with our values for \m.hpi p(j? 

In order to compare with older results in the literature, we have calculated two additional terms, 
both proportional to 7 g, in which the CP-violating weak transition occurs in the meson sector. The 
tree-level kaon-pole contribution to the P-waves will be shown in one of our tables because this is in 
fact the dominant contribution to the commonly quoted result of Donoghue, He, and Pakvasa [6], 
as we discuss below. The one-loop nonanalytic contribution proportional to 7 g occurs at order 
in the chiral expansion. It is related to the model employed by Iqbal and Miller in Ref. [9], and we 
include it here to comment on that result. 

For our numerical calculations, we use the leading-order (in QCD) Wilson coefficients at /i = 

= 1.3 GeV listed in Table XIX of Ref. [30]. In particular, 

y, = -0.096 , (30) 

corresponding to A^ = 325 MeV. This is one of the middle values of in this table, which vary 
from —0.063 to —0.120, depending on the value of A^ and on the renormalization scheme. In the 
rest of this section, we numerically evaluate the weak phases in the A and decays, relegating the 
corresponding evaluation for the S decays to Appendix C. 

The nonfactorizable contributions from Qq to the weak parameters are given by the bag-model 
results in Eq. (25). The resulting s and p amplitudes are collected in Table H, divided by the 
experimental amplitudes of Table I. For the factorizable contributions, the parameters are given 
in Eq. (29) and the corresponding amplitudes are listed in Table HI. In calculating the imaginary 
parts in these tables, we employ the Pq value in Eq. (30), as well as the strong couplings D = 0.8, 
F = 0.5, \C\ = 1.7, and Ti = —2.4. The loop contributions are computed using the results of 
Ref. [35] at a renormalization scale of 1 GeV with fp = f^, and serve as an error estimate of the 
prediction given by the tree contributions. In Table HI, we have separated out the terms containing 


^ This prescription of taking the leading nonanalytic contributions as the uncertainty in the lowest-order amplitudes 
works remarkably well for the real part of the amplitudes. To show this, we use the weak parameters determined 
from htting simultaneously the S-wave amplitudes in Eq. (20a) and the leading-order P-wave amplitudes for 
n ^ AK,Ett provided by Ref. [29] to the measured amplitudes. Thus, hp = 0.49, hp = 1.18, and hp = 1.15, 
all in units of f^Gp . Writing the resulting amplitudes as treeiloop, and excluding 7 g terms, we have 
= 1.25 ±2.28, Ss-^A.- = -1.65 ±2.96, = 0.49 ± 0.92, and = -0.16 ±2.21, all in 

units of Gpm^+. Clearly the corresponding data in Table I are well within these ranges. 
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TABLE II: Ratios of the imaginary part of the theoretical value to the experimental value, for S- and 
P-wave amplitudes, with the weak couplings from contribution only, estimated in the bag model. The 
ratios are in units of rjX^A^. 


Decay mode 

Im Stree 

Sloop 

ImPtree 

ImPioop 

■^expt 

■^expt 

^expt 

^^expt 

X^P'K~ 

-0.09 

-0.09 

-0.15 

0.28 

—!■ Att” 

-0.06 

-0.06 

0.14 

-0.26 


TABLE III: Ratios of the imaginary part of the theoretical value to the experimental value, for S- and 
P-wave amplitudes, with the weak couplings from Qq contribution only, estimated in factorization. The 
ratios are in units of rjX^A?. 


Decay mode 

Im St,ee 

^expt 

Sloop 

Sexpt 

Sloop 

Sexpt 

ImPtree 

Pexpt 

^^expt 

Impoop 

Pexpt 

^^expt 

[I] 

1 > 
i i 

1 

1.13 

1.00 

1.05 

1.10 

0.35 

0.56 

1.33 

-0.66 

0.04 

-0.02 

-0.27 

0.59 

0.61 

-0.28 


Im 7 g. In the P-waves, the 7 g contributions also occur at next-to-leading tree-level order, arising 
from the kaon-pole diagram in Fig. 2. 

In Table IV, we combine the weak phases from the preceding two tables, keeping only the leading- 
order and loop contributions (excluding 7 g terms). We also show in this table another error estimate, 
5(j), obtained from the leading-order amplitudes, but allowing the parameters to vary between their 
tree-level and one-loop values. In making this estimate, we use only the factorization amplitudes, 
as they are are much larger than the bag-model contributions, as seen in the previous two tables. 
Thus, for the P-wave amplitudes, we need the one-loop values of the parameters 6^ p. Employing 
the one-loop formulas for baryon masses derived in Ref. [39], we find 

bp = -0.636 , bp = -0.192 . (31) 


For the P-waves, we note that the factorization parameters in Eq. (26) and the tree-level mass 
formulas 




— ‘^ipD ^f) 




m^-rripf = -|(&z) + 3&^) 




(32) 


— rrij. = —2 (bp + bp) = | (bp — 3bp) , 

derived from Eq. (11), lead to simplihed expressions for the leading-order amplitudes arising from 
the Qe contribution, namely. 


a 


(P) 

Ap7r~ 


GpX4C,fB, 

a/ 2 A/6mg 


(-D - 3P) 




a 


iP) 

E+n7r+ 


0 


a 


(P) 

S“n 7 r“ 



GpX4C,pB, 

a/2 pQrrig 


(-P + 3P) , 


GpXAG.pB, ^^ 

m. 


(33) 
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TABLE IV: Weak S- and P-wave phases from Qq contribution alone, in units of rjX^A^. 


Decay mode 

1 (tree) 

,(loop) 

Vs 

r ,(tree) 

/ (tree) 
(pp 

^aoop) 

c- , (tree) 
0(pp 

A^pTr~ 

1.04 

0.96 

-0.83 

1.18 

0.01 

-0.30 

—> Att” 

0.94 

1.04 

-1.04 

-0.52 

0.33 

0.27 


where the S-decay amplitudes have been included to be used in Appendix C. Consequently, we only 
need the one-loop values of D and F. A one-loop £t to the semileptonic hyperon decays yields [25] 

D = 0.61 , F = 0.40 . (34) 

Using these results, together with their tree-level counterparts in Eqs. (8) and (28), we write the 
ranges 


-0.64 <bj^< +0.18 , 
0.61 < D < 0.80 , 


-0.56 <hF< -0.19 , 
0.40 < F < 0.50 . 


(35) 


We take 5(j)g p to be the largest deviation from 0^^ (in factorization) allowed by these ranges. 

From the nnmbers in Table IV, we may conclude that the uncertainties of (j)g and (jyp are of order 
100% and 50%, respectively, for both decays. This is reflected in onr prediction for the phases, which 
are collected in Table V along with the resnlting phase differences. The errors for the differences 
have been obtained simply by adding the individual errors. We have also collected strong-phase 
differences in the table, from the numbers given in the Introdnction. The errors we qnote in this 
table are obviously not Gaussian. They simply indicate the allowed ranges within our prescription 
to calculate the phases. 


TABLE V: Weak phases in units of rjX^A^, and strong-phase differences, 5g — 6p. 


Decay mode 

4>s 

(j)p 

ds — 4'p 

6g — 5p 

A—fpTr~ 

1.0+ 1.0 

1.2 ±0.6 

-0.2 ± 1.6 

7° ±2° 

—>-A'k~ 

0.9 ±0.9 

-0.5 ±0.3 

1.4± 1.2 

1.1° ±2.8° 


Putting together these results, we finally obtain 


leading to 


A(A°) = Aa = (0.03 + 0.25) A^X^r] , 
A(S:) = As = (-0.05 + 0.13) A^X^ri , 

= Aa + As = (-0.02 + 0.38) A^X^r] . 


(36) 


(37) 
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With the CKM parameter values given in Eq. (17), we have ~ 1.26 x 10 and therefore 


-3 X 10“® < 71^ < 4 X 10“® , 


-2 X 10“^ < < 1 X 10“® , 


(38) 


-5 X 10“^ < 71=^ < 5 X 10“® . 


(39) 


V. DISCUSSION 


We start by comparing our results to those that can be found in the literature. The result most 
frequently quoted is that of Donoghue, He, and Pakvasa [6] given in their Table II, 

yl(A°) = -5 X 10“^ , H(H:) = -7 x lO"® . (40) 

This result was computed using the matrix elements obtained by Donoghue, Golowich, Ponce, and 
Holstein [40]. Recast in the language of our previous sections. Ref. [40] estimated p and 

Im 7 g as the sum of direct and factorizable contributions in the same way we have done in this 
paper. The direct (nonfactorizable) contributions were calculated in the MIT bag model, and we 
agree with their results up to numerical inputs. The factorizable contributions in Ref. [40] are the 
ones they attribute to the quantity We disagree with the calculation of these factorizable 

terms in Ref. [40] in several important ways. 

• For the S'-waves, we obtain a factorizable contribution to hp approximately 4 times larger than 
that of Ref. [40]. This can be traced mainly to a difference in two factors. First, for the chiral 
condensate we use (0|(id|0) = —— —0.012 GeV^, instead of (0|(id|0) = —0.007 GeV^ 
used in Ref. [40]. Second, we employ the value bpiri^ ~ —95MeV in Eq. (27), obtained 
from a hrst-order fit to the baryon-octet masses with Eq. (11), whereas Ref. [40] calculate a 
baryon overlap in the MIT bag model that is equivalent to using bpm^ ~ —43 MeV, with 
rrig = 170 MeV. 

• A second difference in the S'-wave phases (less important numerically) occurs because we use 
hp — O.Shp, as can be seen from Eqs. (26) and (28), whereas the results of Ref. [40] used 
in Ref. [6] correspond to hp = 0. 

• Our most important difference occurs in the P-waves. Our factorization results from leading- 
order yPT calculations arise from the baryon poles. In contrast, the results of Ref. [40] for 
the baryon poles appear to include only the nonfactorizable contributions, and their P-waves 
are instead dominated by the kaon pole, as in Fig. 2. This kaon pole is not included in our 
calculation because it occurs at next-to-leading order in yPT and, moreover, it is further 
suppressed by a factor of rn^/m\ because the pion (and not the kaon) is on-shell. 

We have calculated this kaon-pole contribution (although we do not include it in our hnal 
results) and present it in the sixth column of our Table HI under the heading “Im ■ It 
can be seen from this table that the kaon pole is indeed negligible compared to the baryon 
poles. Studying the calculation of Ref. [40], we believe that their large result for the kaon 


13 



pole is incorrect. The specific error arises in the evaluation of the kaon-pion weak transition 
in the bag model. We show some details in the last part of Appendix A. It is useful to cast 
this issue in the language adopted by the e' literature [30], 


(« 6 > 


0 



m 


1 2 


K 




{Ik 



(41) 


where {Qq)o = (tttt, J = Q\Qq\K). In our estimate, we use a yg corresponding to the value 
= 1 from factorization. For comparison, current lattice estimates are in the range 
= 1 ± 1 [41], whereas the calculation of Ref. [40] is equivalent to ~ 35. 

Despite this disagreement, the numerical value for the P-wave phases based on the results of 
Ref. [40] is similar to ours. This agreement is fortuitous and occurs because the factorizable 
contribution to the baryon poles is roughly equal to 35 times the kaon pole, as can be seen in 
Table III. 


In view of the above, the resulting numerical differences occur mostly in the F-wave phases, 
ours being larger than those found in Ref. [6]. This in turn impacts mainly the phase difference in 
the A case, as 0^p now tend to cancel each other. In contrast, the corresponding phase difference 
calculated using the results of Ref. [40] is much larger (by a factor of 5), being dominated by the 
P-wave phase. In the S case, the two weak phases have opposite signs, and so their difference is 
not suppressed, but instead it is now enhanced (by a factor of 3) with respect to that based on 
Ref. [40]. All these differences lead to the central values in Eq. (36), in comparison to the results 
of Ref. [6] in Eq. (40). An additional problem with the numbers in Eq. (40) is that they follow 
from outdated numerical input for the CKM matrix elements (and also from the use of the large 
old value (5| ~ —18° for the S decay [42]). 

Next we turn our attention to the vacuum-saturation calculation of Ref. [7]. Our results in 
Tables II and III indicate that the factorization contribution is significantly larger than the direct 
contribution to the S- and P-wave phases. For this reason, we would expect our results to agree 
with those of Ref. [7] in which the direct contributions are ignored. We find that we agree with the 
value of the S'-wave phases up to numerical input, but that we disagree with the value of the P-wave 
phases in Ref. [7]. This disagreement is easy to understand. Our P-wave phases are dominated 
by the baryon-pole contribution, whereas in Ref. [7] only the kaon-pole contribution is included. 
The vacuum-saturation calculation of the kaon pole, corresponding to = 1, is a significant 

underestimate for the P-wave phases as seen in Table III, where the kaon pole corresponds to the 
column labeled “Im p|ree”. 

To summarize then, the bag-model calculation of Ref. [40] significantly overestimates the contri¬ 
bution of the kaon pole to the P-waves and apparently misses the important factorization contri¬ 
bution of the baryon poles, although accidentally results in P-wave phases numerically similar to 
ours. Furthermore, it underestimates the A-waves and therefore yields an asymmetry dominated by 
the P-wave phase. The vacuum-saturation calculation of Ref. [7] misses the dominant baryon-pole 
contribution to the P-wave phases and results in an asymmetry dominated by the phase of the 
S'-wave. In our complete calculation at leading order in yPT, the phases of the S- and P-waves 
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are comparable, and in the A case this leads to a smaller central value for the predicted asymmetry 
(the two phases tend to cancel). 

It is difficult to place the calculation of Ref. [9] in our framework due to signihcant technical 
differences in the evaluation of loop integrals. Nevertheless, there is a rough correspondence between 
that calculation for the S'-waves and the terms in Table III labeled “Im . In our final results, 
such terms appear in the quoted uncertainty because they are part of the subleading amplitudes 
that cannot be calculated completely at present. 

In conclusion, we have presented a complete calculation of the weak phases in nonleptonic hy- 
peron decay at leading order in heavy-baryon chiral perturbation theory. We have estimated the 
uncertainty in our calculation by computing the leading nonanalytic corrections. We have compared 
our results with those in the literature, pointing out several errors in previous calculations. To im¬ 
prove upon the results presented in this paper, it will be necessary to have a better understanding 
of the P-waves in nonleptonic hyperon decay. 
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APPENDIX A: BAG-MODEL PARAMETERS 


In this appendix, we summarize the derivation of the formulas in Eq. (24), which describe the 
nonfactorizable contributions to the weak parameters pQ, estimated in the MIT bag model.^ We 
also provide the numerical values of the parameters a and b in these formulas. Lastly, we evaluate 
the kaon-pion matrix element of the leading penguin operator in the bag model. 

Assuming a valence-quark model of baryons, using the totally antisymmetric nature of their color 
wavefunctions and the relations [30] 


Qa Ql T Q 2 ffi Qs > Q 9 2 2 ^3 ’ 

one hnds for baryons B and B' 

{B'\Q,\B) = -{B'\Q,\B) = {B'\Q,\B) = , 

(B'lOs + Qt\B) = {B'\Q, + Q,\B) = {B'\Qj + Og|B) = 0 . 


(Al) 


(A2) 


^ An introductory treatment of the bag model can be found in Ref. [22] 
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Therefore, only need to be evalnated. For the parity-conserving parts of < 5 i^ 5 ^ 7 , we 

derive the bag-model matrix elements^ 

(n|Qi|A) = (nlQslA) = - 2 (n|Q 7 |A) = -^ 6(0 + 6 ) , 

(A|Q,|H»)=2v^(a + ()), (A|Q,|H») =-2(A|Q,|H“) = , (A3) 

= 0 , = - 2 (E-|(J,|fi-) = -4V3(a - it) , 

up to factors of where a and b will be described shortly. From these results and Eq. (23), 

we then obtain 

Ad + SAd = ^^(c, - C, + C, - C, + C, - C,„ + C, - Q - iC, + iC*) 6{-a - b) , (A4) 

Ad - 3Ad = [(C, - C, + C3 - C. + C, - C,„) 12(0 + A) + (Cj - C5 - iC, + iCs) 16 a] , 

(A5) 



^ In keeping with Eq. (23), we have excluded from these results the 27-plet components of y and the (SljSj^) 
component of Qy, the strong penguin operator being purely (Sl, Ij^). Furthermore, in the 0-S* matrix-elements 
we have taken into account the fact that the spinors for decuplet baryons in the chiral Lagrangian are spacelike [23], 

■Uh. • Ujy < 0. 
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where 


M = 


p 


a/ {2uj^ — 2uj + mR) sin^p 


p = \/up' — im?R^ , 


e = 


LO 


mR 


uj + mR ’ 


(AlO) 


with ui being determined from tanp = p/{l — uo — mR) and m the quark mass in the bag. Numer¬ 
ically, following Refs. [40, 43], we take R = 5.0GeV“^ for octet baryons and R = 5.4GeV“^ for 
decuplet baryons. Since the weak parameters belong to a Lagrangian which respects SU(3) 

symmetry [£^ in Eq. (13)], in writing Eqs. (A3) and (A7) we have employed SU(3)-symmetric 
kinematics.^ Specihcally, we take m = 0 for all quark flavors. Thus, we hnd for octet baryons 


a = 1.40 X 10-^^GeV^ , 

b = 0.64 X 10-3 GeV3 , 

(All) 

and for decuplet baryons 



a = 1.11 X 10“^GeV^ , 

b = 0.51 X 10-3 GeV3 . 

(A12) 


Finally, we evaluate the K-io-n transition in the bag model, which occurs in the kaon-pole result 
of Ref. [40], as discussed in our Sec. V. The matching of the dominant |AJ| = 1 part of the weak 
Hamiltonian in Eq. (14) to the weak chiral Lagrangian in Eq. (13) involves in this case 

{k-{ p)\{n^)^\K-{p)) = -278• (A13) 

Goncentrating on the Qq contribution alone, we hnd the bag-model matrix element 

{'k-\Q,\K-) = -12(a + 6) (A14) 

where the factor arises from the normalization of the bag states for the mesons [22, 40]. 

From the preceding two equations, we obtain the contribution 

= -^Gg6(a-h6) . (A15) 

To determine the values of a and b in this equation, we use R = 3.3 GeV“\ after Ref. [40], and 
again set m = 0 for all quark havors. It follows that here 

a = 4.87 X 10-=^GeV=^ , b = 2.23 x 10“^ GeV^ . (A16) 

At this stage (in their equivalent calculation). Ref. [40] proceeds by setting = m^ and AR^R^ = 
2mj^. As a consequence, 

Im 78-^0 = 6(a + &) ^^°^-^ = 638 pg (A17) 

v2 

in units of \/2 f^G pm^+ with the Bq value in Eq. (28). Gomparing this result with Eq. (29) 

then indicates that the bag-model calculation of Ref. [40] yields ~ 35, which is unacceptably 

large. 


^ We note that in the SU(3)-symmetric limit the bag parameters above are related to the parameters A and B of 
Ref. [40] by H = 47r(o - b)R^/M^ and B = SirbR^/Af'^. 



17 



APPENDIX B: WEAK PARAMETERS IN FACTORIZATION 


To derive the factorizable contributions to the imaginary part of the parameters pQ, we start 
from the observation that the quark-mass terms in the QCD Lagrangian can be written as 


2Br 




(Bl) 


where — 75 )^ and q^= |(1 + l 5 )q, with q = {u d s)"’". It follows that 


■Q'iL 5'fcR ~ 2 -Bq 




Sx. 


Ik 


'%R^kL — 2-^0 I 

^X'lk 


(B2) 


Then, using in Eqs. (11) and (12), we have the correspondences 


-QiLQkR ^ + + bp{^^BB^^ + ai:l{BB) 

+ + 5/%d, 

+ aBo-tqS'-E'^EEt)^, + ... , (B3) 


-%RlkL ^ + + aB„{BB) 

+ 2B„L,(Ea'‘Eta,.E)^, + ... , (B4) 

where the ellipses denote additional terms from CA that do not affect our result. Consequently, 
for the penguin operator 


we obtain 


Qe = Y 1 ^(1 + 75)99(1-75)5 = ~^ Y 1 ^ l 9 r 9 r 5 l , 

q=u,d,s q=u,d,s 


-|«8 ^ fB„(bo(B{(^h(,B}) + b^{B [«th«,B]» 

+ PB„cf‘ehiT„ + 2pBf,L,(h8^E8>‘B') + 


(B5) 


(B 6 ) 


where only the terms that correspond to leading-order chiral perturbation theory have been shown. 
Comparing this expression with the weak Lagrangian in Eq. (13), we then infer that the contribu¬ 
tions of a factorized Qq to the weak parameters are 


k _ f2 

rip — 


V 2 


hc= 




8 Q rB, bp , 

8C,fB,c, 




72 


8C,rB,bp, 


78 = ^IQC,B^,L,. 


(B7) 


72 
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APPENDIX C: CP-VIOLATING ASYMMETRIES IN S ^ ATtt DECAYS 


The S'-wave amplitudes in S —Nn can be expressed in terms of their components 2 /’ 

where the I in the second subscript denotes the isospin of the Ntt state. Thus we have® 


’S+^n7r+ 


i ( 2 s,^ e‘'^ii + ^31 e'^3i) e'^? + i e‘^13 


_ Q pl'Ass I pJU3 

5 ^33 ^ J ^ > 


Sr*-,,. = 5^(2S„e'*?.+S3ie«)e“f-f (^S,3e‘<’&-2^S33e'*?.je“3 , (Cl) 




(s,3e'*& + \/iS33e‘*t.) e“l , 


where S 21 and (1>2\ai\ 21 strong iVTT-scattering and weak CP-violating phases, respectively, 

and |A/| = | components have been ignored. The P-wave amplitudes can be similarly expressed. 
For each of these decays, one can construct the counterpart of the CP-violating asymmetries 
using [6] 


T a + T a 
T a — f d 


(C2) 


One then has 


.4(S+) ^ A 


S+^n7r+ 




sin((5f - 5f) sin(0f “ ) + ^ sin((5f - (5f) sin(0f - 0f) 


P, 




+ ^ sin((5f - (5f) sin(0f - 0f) + sin((5f - (5f) sin(0f - 0f) 

cos((5f - (5f) + ^ cos((5f - <^ 3 ) + ^ cos((5f - 5f) + cos(5f - (5f) 


, (C3) 


® In the phase convention that we have adopted to write down these amplitudes, the isospin states Idjp) for the 
hadrons involved are |i;+) = —11,1), |I]“) = |1,—1), |p) = |l/2,l/2), \n) = |l/2,—1/2), |7r+) = —11,1), 
|7r°) = |1,0), and |7r“) = |1,—1), which are consistent with the structure of the (f and matrices in the chiral 
Lagrangian. 
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^(^o) — ^S+^PttO 


sin((5f — (5f) sin( 0 f 

P 3 


sin 


COS 


iin(5f - <^f) sin( 0 ^ - 0 f) + sin(5f - (5f) sin( 0 ; 

S 3 P, 


Pi) 


Pi 

^ cos( 53^ - 5f) + cos( 53^ - 5f) , (C4) 


- '^f) -1 - <^ 3 ") - ^ 


A(S:) ^ 


- tan((5| 


■<^3j 


where 


= Sn+^^S., 


31 > 


iS — ‘S'll 011 
^ ■ ^1 


1 _|_ ^31 _|_ - ~^33 _|_ ^ P33P33 

^V5^33^V5 A3^5^,3A3 


^3 = ^13 


2 C 

5 ^33 > 


iS — Pl3 013 
^ 


(C 6 ) 


the P-wave counterparts being similarly dehned, and the weak | A/| = | phases have been neglected. 

To estimate the weak phases, we follow the prescription proposed earlier, obtaining the real 
part of the amplitudes from the values extracted from experiment under the assumption of no CP- 
violation and calculating the imaginary part from the leading-order amplitudes in Eq. (20) with 
the values of Imhj^ p provided in Section III. To hnd the real part, ignoring the strong and weak 
phases, we hrst derive from Eq. (Cl) 


^3 = s 


s+- 


ir+ ^^2 "^S+^pTfO ) 


‘S '33 “ 


18 I '®S+^n7r+ 


3 S '33 — Ss+^r^ 7 r^ 

0 — 


V 2 + 2 S 5 


(C7) 


and analogous expressions for the P-waves. From the experimental values in Table I, we then 
extract, in nnits of Gpm^+, 


S-^ = -0.95 ± 0.04 , ^33 = 1.95 ± 0.02 , ^33 = -0.11 ± 0.04 , 

(C 8 ) 

Pi = 2.64 ± 0.04 , Pi 3 = 0.01 ± 0.03 , P 33 = -0.11 ± 0.05 , 

The imaginary part of the amplitudes are obtained using Eqs. (20), (25), and (29), as well as the 
isospin relation v^Os+^pttO ~ ®s+^n 7 r+ ~ ®s-^n 7 r- l^-^l ~ 3 dominance. Thus, we have in 

units of 


^ = (-0,04+ 1.02)+ (-0.03 + 1.30) , 

IlXl P 

= (0.02+ 0.09)+ (-0.07+ 0.31) 

-i 1 


ImPi3 



(-0.04 + 0.99) + (-0.03 + 1.26) 


ImPi3 


pexpt 

■^13 


(7-41) + (-15-59) , 


(C9) 
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where the numerators on the left-hand sides are the central values in Eq. (C8), and we have written 
each result as (tree)-|-(loop), with the two numbers within each pair of brackets being bag-model 
and factorization contributions, respectively. In Table VI, we collect the weak phases resulting from 
these ratios. 

We also show in Table VI another error estimate, <50, obtained from using the leading-order 
amplitudes and allowing the parameters to vary between their tree-level and one-loop values, as 
discussed in Sec. IV. In making this estimate, we again employ only the factorization contributions 
[for the P-waves, we use the S amplitudes in Eq. (33)], which are are much larger than the bag-model 
ones, as seen in Eq. (C9). 


TABLE VI: Weak S- and P-wave phases in S ^ Ntt decays from Qq contribution alone, in units of riX^A^. 


j S.tree 

4>i 

^S.loop 


i S.tree 
^13 


c iS.tree 

0013 

j P.tree 

(Pii 

iPAoou 

c / P.tree 

Hii 

i P.tree 

^13 

j P.loop 
<l>13 

r- (P.tree 
0013 

0.98 

1.27 

-1.65 

0.95 

1.23 

-1.61 

0.11 

0.24 

-0.05 

-34 

-74 
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We may, therefore, conclude that the uncertainties of the weak phases are all of order 200%. This 
is reflected in our prediction of the phases, which are collected in Table VII. The corresponding 
strong phases have been measured [11] and their values have also been included in this table. 


TABLE VII: Predicted weak phases, in units of and measured strong phases. 



<^f3 


dfa 


<51 

sf 

8^ 

1.0 ±2.0 

1.0 ±2.0 

0.1 ±0.2 

-40 ± 80 

9.4° ± 1.0° 

-10.1° ± 1.0° 

-1.8° ± 1.0° 

-3.5° ± 1.0° 


From the central values of the isospin amplitudes and the phases in Eq. (C8) and Table VII, 
respectively, we obtain 

A(S+) = 3.9 X 10“^ , A(S+) = 3.6 x 10“® , A(S:) = -8.3 x 10“^ , (CIO) 

where we have used A^X^rj = 1.26 x 10“^ as before. In this case our estimate is a very rough 
one, as its uncertainty is larger than those for the other hyperons. This is due to the (apparently 
accidental) smallness of P^g and its large experimental error, indicated in Eq. (C8), as well as to the 
already sizable uncertainties quoted in Table VII. In order to have a more quantitative estimate of 
the uncertainties, these modes will have to be revisited when better measurements of the amplitudes 
become available. 
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